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ABSTRACT 

Galactic bulges are known to harbour central black holes whose mass is tightly cor¬ 
related with the stellar mass and velocity dispersion of the bulge. In a hierarchical 
universe, galaxies are built up through successive mergers of subgalactic units, a pro¬ 
cess that is accompanied by the amalgamation of bulges and the likely coalescence 
of galactocentric black holes. In these mergers, the beaming of gravitational radiation 
during the plunge phase of the black hole collision can impart a linear momentum kick 
or “gravitational recoil” to the remnant. If large enough, this kick will eject the rem¬ 
nant from the galaxy entirely and populate intergalactic space with wandering black 
holes. Using a semi-analytic model of galaxy formation, we investigate the effect of 
black hole ejections on the scatter of the relation between black hole and bulge mass. 
We find that while not the dominant source of the measured scatter, they do provide 
a significant contribution and may be used to set a constraint, Ukick ^ 500kms“^, 
on the typical kick velocity, in agreement with values found from general relativistic 
calculations. Even for the more modest kick velocities implied by these calculations, 
we find that a substantial number of central black holes are ejected from the progen¬ 
itors of present day galaxies, giving rise to a population of wandering intrahalo and 
intergalactic black holes whose distribution we investigate in high-resolution N-body 
simulations of Milk-Way mass halos. We find that intergalactic black holes make up 
only 2 — 3% of the total galactic black hole mass but, within a halo, wandering 
black holes can contribute up to about half of the total black hole mass orbiting the 
central galaxy. Intrahalo black holes offer a natural explanation for the compact X-ray 
sources often seen near the centres of galaxies and for the hyperluminous non-central 
X-ray source in M82. 

1 INTRODUCTION 

The discovery that galactic bulges harbour supermassive 
black holes whose masses are correlated with the properties 
of the bulge suggests a close connection between the forma¬ 
tion of galaxies and the formation of black holes. The mass 
of the galactocentric black hole vari es approximately linearly 
with the stellar mass of the buf^ejKormen^^iRichstony 
IIQQItI : iMagorrian et all Il998l : lMcLure^^DunioD[]200^ and 
ronghly as the fourth power of the bulge velocity dis¬ 
persion j jFen^esefc MerrittI l200fi . lOebhardt et al.l 120001 
lTremain^^^T l200^~ Correlations wit h the near infrared 
bulge luminosityjMarcon^&^^n^^OO^Iand with the bulge 
light concentration llOrEdian^^^ J|2QmD have also been 
found. These correlations are remarkable because they link 
phenomena on widely different scales - from the parsec scale 
of the black hole’s sphere of influence to the kiloparsec scale 
of bulges - and thus point to a connection between the 
physics of bulge form ation and the physics of blac k hole 
accretion and growth llMilosavlievic fc Merritt|l2001^ . The 
simplest interpretation is that both, black hole and bulge 


growth, are driven by the same process whose nature, how¬ 
ever, remains unclear. 


Various models for the growth of black 

holes in galaxies have been studied (e.g. 


Haehnelt & R.eesI Il9^ 

Silk & R.eefJ 1199(1: 

Cattaneo. Haehnelt Kr, R.eesJ Il999l: 

Kauffmann & HaehnehJ 

200(1 lOstrikerl 1200(1: IVolonteri. Haardt & MadaiJ l2002l: 

Di Matteo. Soringel & Hernouistl l200!l and others). In 


the context of a hierarchical cold dark matter universe, 
a plausible explanation for the tight correlation be¬ 
tween bulge and black hole properties is that the galaxy 
mergers or disc instabilities that induce bulge growth 
via b ursts of star formation also f ee d the central black 
hule_j [KaA^mja^ fc Haehneld 1200(1 : ICroton et alJ l200fil : 
iBower e^n]20oT ' A simple implementation of this model 
follows from assuming that, as cold gas condenses into stars, 
a certain percentage of the gas is forced into the centre of 
the galaxy and accreted by the black hole. Models based 
on this and related prescriptions successfully reproduce the 
rriBH — rrih relation (hereafter we use this term to refer to the 
relation between black hole mass and stellar bulge mass) 
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JCroton et alJ l200fil : iBower et alJ I200RI : iMalbon in pred 
1200,'^) . 

An interesting aspect of the correlations between black 
hole mass and bulge properties is that they seem to apply 
over a range of hve to six orders of magnitude in black h ole 
mass (iTremaine et al ] l2002HGebhardt. Rich fc Holl2002^ . If 
a simple model of the kind just mentioned for the simul¬ 
taneous growth of black holes and bulges is correct, then 
there are two direct consequences which we explore in this 
paper. The hrst is that black holes should exist in bulges 
of all luminosities including dwarf ellipticals and satellites 
of brighter galaxies like the Milky Way. The second is that 
black holes will likely merge as their hosts bulges collide. 

Binary blac k holes orbiting each ot her emit grav¬ 
itational waves jPeters fc Matthe^ Il96.?ll . Using quasi- 
Newtonian methods to study the orbita l decay due to grav¬ 
itational wave emission. lFitchet3 Jl98,^l found that the sys¬ 
tem will eventually enter a plunge phase, causing the black 
holes to coa lesce emitting a burst of gravitational waves. 
IPeresI (ll962^ found that, in addition to transferring energy 
out of the emitting system, gravitational radiation can also 
take with it linear momentum. As a result, the centre of 
mass of the system recoils in a direction specihed by the 
boundary conditions of the last stable orbit. 

The astrophysical implication of this linear momen¬ 
tum kick (the “rocket effect” or “gravitational recoil”) 
is that black holes may be ejected from galactic bulges 
if the potential is shallow and the kick is large enough 
jMndau fc Quataertll20oS iMerritt et alJ 1^041 : lEnoki et a.lJ 
l20n4l . In theory, this could lead to a sizable population of 
extragalactic black holes which could, in principle, dominate 
the black hole mass function. Our aim in this work is to 
examine the importance of such kicks for the galactic black 
hole population. In particular, we consider the scatter on the 
msH — mb relation in an attempt to constrain the relatively 
uncertain kick velocity, as well as the nature and spatial 
distribution of a possible extragalactic population of ejected 
black holes. We model the growth of bl ack holes us i ng th e 
semi-analytic galaxy formation model of ICole et al.l ll200(il . 
using two methods for obtaining merger trees: Monte Carlo 
techniques and high resolution N-body simulations (in which 
the trajectories of recoiling black holes can be tracked.) 

This paper is organised as follows. In section |21 we re¬ 
view the physics of the gravitational recoil. In section we 
describe how we model the formation and ejection of black 
holes in our semi-analytic model. In section |3] we deter¬ 
mine the effect of black holes ejected from the progenitors 
of present day galaxies on the mnn — ms relation. In sec¬ 
tion 1^ we use a set of high resolution N-body simulations 
of galactic halos to track the location of ejected black holes. 
We conclude in Section|^where we discuss the possible con¬ 
sequences of an extragalactic black hole population. 


2 THE PHYSICS OF KICKS 

The emission of gravitational radiation is a gen eric feature o f 
any massive asymmetrically collapsing system llPere iim. 
In order to be accompanied by a linear mom entum kick, th e 
gravitational radiation must be asymmetric llFitchettll983ll . 
For two coalescing black holes of unequal mass this occurs 
as the gravitational radiation from the lighter, more rapidly 


moving partner is more strongly beamed. As two black holes 
orbit each other, gravitational radiation causes their orbits 
to shrink and circularise until the last stable circular or¬ 
bit is reached. Using perturba tion theory in the weak field 
approximation, iFitchetil (Il98.8h calculated the kick velocity 
from this final orbit to be: 

= 1480 ( T«l) ( MajhfalTy k„,-. (1) 

where risco is the radius of the innermost stable circular 
orbit, q = mi/m 2 is the mass ratio with the convention 
m 2 >mi,and/(q) = g^(l — 5 )/(l-|-(?)® which reaches a max¬ 
imum, /max = 0.0179, for a mass ratio of qmax ~ 0.382. Un¬ 
fortunately, the weak field approximation used by Fitchett 
becomes invalid as the binary approaches the plunge phase 
and it is here that the contribution to the recoil velocity 
becomes largest. 

Two decades later, Favata, Hughes & Holz (2004) 
used perturbation theory to show that, in the presence of 
strong fields, gravitational kicks ar e not expected to exc eed 
OOOkrns”^. In a companion paper. iMerritt et alJ (|20o3) in¬ 
cluded the effect of the larger partner’s spin and obtained 
an upper limit of ^ 500kms“^. More recently Blanche!, Qu- 
sailah & Will (2005) performed a higher order calculation of 
the recoil for the special case of the coalescence of two non¬ 
rotating black holes, obtaining an upper limit of 300 kms“^. 
However, none these calculations apply to all mass and spin 
ratios and kick velocities as large as 1000 km s“^ cannot be 
definitively ruled out. 

In light of the uncertainty in the recoil velocity, we sim¬ 
ply assume that the kick is directly proportional to Fitchett’s 
scaling function f{q) and write the recoil velocity in terms 
of a prefactor velocity Upf, 

Ukick = Vpt . (2) 

/max 

We then allow Vpi to vary between 0 kms“^ (i.e. no kick) 
to 1000 kms”*^ with the aim of constraining this parameter 
empirically. 


3 THE GALFORM SEMI-ANALYTIC MODEL 

In this section, we briefly explain how we use the semi- 
an alytic galaxy form a tion model GALFORM descr ibed iu detail 
by j^le et al.l (I 2 OOOII . iBenson et al.l ^2^^2^ and iBower et alJ 
^2005^ to model the growth and evolution of galaxies. The 
growth of dark matter halos by mergers is encoded in a 
merger tree. Along each branch of the tree, the following 
physical processes are calculated: (i) shock-heating and viri- 
alization of gas within the gravitational potential well of 
each halo; (ii) radiative cooling of gas onto a galactic disc; 
(iii) the formation of stars from the cooled gas; (iv) the ef¬ 
fects of photoionization on the thermal state and cooling 
properties of the intergalactic medium; (v) reheating and 
expulsion of cooled gas through feedback processes asso- 
ciated with stellar winds and supernovae explosions (see 
iBenson et alJliooifi : (vi) the evolution of the stellar pop¬ 
ulations; (vii) the effects of dust absorption and radiation; 
(viii) the chemical evolution of the stars and gas; (ix) galaxy 
mergers (which, depending on the violence of the merger, 
may be accompanied by starbursts and the formation of a 
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bulge - see lBaugh et, alJll998h : (x) the evolution of the size 
of the disc and bulge. 

We used two different methods for obtaining dark 
matter halo merger trees. In the first instance, we con¬ 
structed the trees using a Monte-Carlo algori t hm to de¬ 
term ine merger rates (e.g. iLacev fc Cold ll99A ICole et ^ 
I 2 OOOII . In the second instance, we explicitly extracted the 
merger trees by followi ng the merging o f dark matter halos 
in N-body simulations (iHellv et all20f).‘fl . The latter method 
allows us to apply the semi-analytic formalism directly to 
the simulations. More details are given in Section 10 


3.1 Discs 


Stellar discs are formed as rotating diffuse halo gas cools and 
settles in a halo. In the GALFORH model, discs are assumed 
to have an exponential surface density profile, 

E(r) = Enexp , (3) 

where ru is a characteristic disc length and Ed is the central 
surface density. The potential of such a mass distribution as 
a function of just radius (i.e. in the plane of the disc) is 


= 0) = -7vGEor[Io{y)Ki{y) - Ii{y)Ko{y)] (4) 


where y = rj^TY). The functions In(y) ^-nly) are mod- 
ified Bessel functions of the first and second kind (e.g. 
iBinnev fc Tremainelll98^ . Differencing the potential at r = 
0 and r = oo, we find that the escape velocity from the 
centre of the disc is 


2 

^D,esc 


3.36 


GMo 
r'D,i/2 ’ 


( 5 ) 


where >" 0 , 1/2 is the half mass radius of the disc. 


3.2 Bulges 

In the GALFORH model, the principal mechanism for forming 
both elliptical galaxies and the bulges of luminous spiral 
galaxies is galaxy mergers. When two galaxies of comparable 
mass merge, this is termed a major merger and all the stars 
from both merging partners are assumed to form a single 
spheroid while any gas present is consumed in a burst of 
star formation. In minor mergers there is no burst of star 
formation; only the stars from the smaller galaxy are added 
to the bulge of the larger galaxy, while the gas is added to 
the disc of the larger galaxy. 

In the context of the GALFORH model, ICole et ll200ff) 
also considered the possibility that dynamical instability in 
selfgravitating disc s could result in the form ation of galac¬ 
tic bulges (see also IMo. Mao fc Whft^ lIlOOjjB . They found 
that this mechanism was unlikely to contribute significantly 
to the formation of large bulges, but could be important in 
creating the bulges of lower luminosity galaxies. The crite¬ 
rion they used to determine whether a cold disc was unstable 
was = 1-1, where 

_ Vjnax 

“ (GMdisc/rD,i/2)i/2 ^ ^ 

iEfstathiou. Lake fc Nee^roDontell982^ . Here Knax is the cir¬ 
cular velocity at the disc half-mass radius, >" 0 , 1 / 2 . We adopt 
the slightly lower value of = 1.05 in order to obtain 


a distribution of bulge to total luminosities similar to that 
found in the Sloan Digital Sky Survey by iTasca fc Whitel 
l|200^. For such unstable discs it was assumed that all the 
disc stars would be transformed into bulge stars and any 
cold disc gas would undergo a burst of star formation. The 
assumption that bulges, and as we discuss later, black holes 
grow as a result of discs becoming unstable is also assumed 
in the AGN feedbac k models incorporated in th e latest semi- 
analyt ic models of ICroton et al.l ll200,‘^l and iBower et alJ 
yOO^. To illustrate the importance of this mechanism, we 
present models both with and without this secondary route 
for generating galactic bulges. 

The densitv orofile we adoot for bulees is thelHernouisti 
CH) model. 


Mb a 1 
2n r (r -I- a)® ’ 


( 7 ) 


where a is a scale length and Mb is the total bulge mass. 
The potential of such a profile is found by integrating the 
Poisson equation and the escape velocity from the centre of 
such a spheroid is given by 


2 


Vb, 


esc 


4.83 


GMb 

>"B,1/2 ’ 


( 8 ) 


where >"b,i/2 is the bulge’s half mass radius. 


3.3 Black Holes 

We assume that the growth of black holes is proportional 
to the growth of bulges. Whenever an event that adds stel¬ 
lar mass to the bulge occurs, we assume that a fraction, 
/bh, of that mass is accreted by the central black hole. 
This assumption is motiva ted in part by models of blac k 
hole growth such as that of iKauffmann fc HaehnelS J200Ctl . 
wherein black holes accrete a proportion of the cold gas be¬ 
ing consumed in a burst of star formation. In our model, 
black holes grow both as a result of such star formation 
bursts and also by the accretion of stars in mergers or as 
the result of disc instabilities. We assume that in each case 
the same mass fraction, /bh, is channelled onto the black 
hole. This simplifying assumption has the virtue that, in 
the absence of black hole ejection, all bulges will have black 
holes that sit precisely on the tubh — mb relation and so 
we can, in principle, use the scatter induced by black hole 
ejections to set a hrm upper limit on the kick velocities. The 
precise value of /bh ne eded to match the data has changed in 
the li terature (e.g. see iMagorrian et al.1ll998l : iMerritt et alJ 
[ 2 OO 4 I) and we choose fsH = 0.001, the value published 
by iMcLure fc Dunlod feoOSll from estimates of black hole 
masses for 72 active galaxies. 

During a merger, we assume that a fraction of any new 
stellar bulge mass is added to the existing black hole, 

Mbb, new — A^BH, old + /bh AMbulge. (9) 

If the black hole is ejected during the merger, we assume that 
a new black hole is born whose mass is equal to the mass 
fraction that would have been added to the preexisting black 
hole had it not been ejected, that is: 

Mbb, new = /bh AMbulge. (10) 

This black hole then becomes the focus of mass accretion 
during subsequent episodes of bulge growth. In this way. 
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we ensure that all bulges contain a black hole at all times 
and that, in the absence of gravitational recoil, all black 
holes at all times have masses directly proportional to their 
bulges and hence lie on a perfect msH — rrih relation with 
zero scatter. The implicit assumption that black holes exist 
only in bulges and not in discs, and that black holes grow 
during mergers rather than during quiescent star formation 
is motivated simply by the desire to reproduce the empirical 
mBH — rrih relation. 

We assume that once two galaxies have merged, their 
central black holes coalesce instantaneously. This should be 
a good approximation as the black holes will merge on a dy¬ 
namical timescale while the time between galaxy mergers, 
which is determined by the hierarchical growth of structure, 
is typically much longer. Note that even when the merger 
rate of dark matter halos is relatively high, the galaxies or¬ 
biting inside the dark halos will only merge once their orbits 
decay due to energy loss by dyna mical friction against the 
halo material (iLacev fc ColellT99.?ll . This implies that three- 
body mer gers in which ejections could occur v ia a slingshot 
effect (e.g. ISaslaw. Valtonen. fc Aarset3ll974l are expected 
to be rare. 

The mass ratio of the merging black holes determines 
the recoil velocity, Ukick , according to equation (|^. Thus, 
to determine whether the recoiling black hole escapes from 
the galaxy we add the bulge and disc potentials and require 

Wick > Wsc = 3.36-h 4.83-. (11) 

’'D,1/2 ’'B,1/2 

Note that we neglect the contribution of the dark matter 
halo. This is a good approximation over the lengthscale of 
the galaxy because the dark matter potential has a very shal¬ 
low gradient. However, the dark matter potential becomes 
increasingly important for black holes that escape the galaxy 
and, in Section EH we use an alternative method to follow 
their orbits. If the recoil velocity, Ukick, is less then Vesc we 
assume that the black hole will not escape but rapidly return 
to the centre of the galaxy. This should be a good approxi¬ 
matio n: according to the calculations of iMadau fc OuataertI 
l|20^), recoiling black holes that have insufficient kinetic en¬ 
ergy to escape undergo damped oscillations about the galac¬ 
tic nucleus which decay in a few dynamical times. 

Our modelling of black hole growth thus automatically 
reproduces an empirical rriBH — mh relation. If the prefactor 
velocity, Upf, in equation (|2|l is zero, the black hole masses 
are all directly proportional to the bulge mass of their host 
galaxies. If Upf is not zero, gravitational kicks introduce scat¬ 
ter in the rriBH — nib relation which we now investigate and 
compare with observations. 


4 THE EFFECT OF VELOCITY KICKS ON 
THE ttibh — nib RELATION 

In order to investigate how the tubh — nih relation is af¬ 
fected by the ejection of recoiling black hole merger rem¬ 
nants we have studied a sample of 1000 halos of final mass 
IO^^Mq. Merger trees for each halo were generated u sing 
the Monte-Carlo method described bv ICole et alJ (l200(t> and 
the GALFORM rules of galaxy formation were applied to each 
branch of the tree. The growth of each central black hole 
(which tracks the growth of the bulge) was calculated as 





Figure 1. The z = 0 rriBH — relation for the case in which 
black hole merger remnants have the kick velocity given by equa¬ 
tion m with Upf = 300 kms The diagonal line represents the 
‘ideal’ rriBH ~^b relation for which MBn/Mbuige = /bh = 10“®. 


described in Section 3.3. Our fiducial model includes the 
process of bulge formation by disc instability discussed in 
Section 3.2. For every black hole merger we consider kick 
velocities corresponding to values of the prefactor, Upf, in 
eqn (2) in the range 0 - 1000 hms”*^. 

The ejection of a black hole from the bulge is a source 
of scatter in the ttibh — rrib relation. Fig.^shows the 2 = 0 
rriBB — mb relation for 1000 halos assuming a moderate pref¬ 
actor kick velocity of 300 kms“^. Note that black hole mass 
may only be scattered downwards by ejections. Galactic 
bulges that have never experienced an ejection lie exactly 
on the diagonal line. 

In the absence of kicks, the model, by definition, gives a 
perfect ttibh — mb relation with no scatter. As the kick veloc¬ 
ity increases, a tail of bulges that host black holes of reduced 
mass appears, while the remaining galaxies whose progeni¬ 
tors never lost their black hole still lie precisely on the rela¬ 
tion. The form of the scatter about the itibh — mb relation 
and its dependence on the prefactor Upf are shown in Fig.|21 
The upper two panels show how deviations from the ideal 
relation depend on bulge mass for both Upf = 300 kms“^ 
and 500 kms“^, while the lower two panels show histograms 
of these deviations averaged over all bulge masses. We see 
in these lower panels that the distribution of deviations 
from the ideal relation is very non-Gaussian, with a tail 
extending to very low values of logio(MBH/iMideai,BH)- For 
Upf = 300 kms“^, 2.6% of bulges have black holes with 
mass smaller than 20% of the ideal mBH — mb value. For 
Wpf = 500 kms“^ this fraction jumps to 18.9%. The rms 
width of this distribution is not a good statistical descrip¬ 
tion of the scatter in the mBH — mb relation since the rms 
is dominated by the noisy tail. We use instead a more ro¬ 
bust measure of width. We calculate the widths that contain 
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10glo(MB„/M,deal,Bn) V [kniS"'] 


Figure 2. Top: residuals from the ruBH — relation. Plotted 
here is the dependence on bulge mass of the ratio of the black 
hole mass to the value for the ‘ideal’ rriBH ~ relation. Bottom: 
histogram of deviations from the ideal rriBH ~ "^b relation. 


68 % and 86% of the distribution, and then define cres and 
as6 to be half and one third of these widths respectively. For 
a Gaussian distribution both of these measures would equal 
a, the usual rms width. Fig. El shows how these measures 
of scatter vary as a function of Upf and compares them to 
published observational determinations. The non-Gaussian 
nature of the scatter results in the two curves being signifi¬ 
cantly different. 

We immediately see that the contribution to the scatter 
in the rriBH — rnh relation from ejected black holes is very 
sensitive to assumptions concerning disc stability. For the 
case where unstable discs are assumed to form bulges, large 
values of Upf are strongly ruled out as they would produce 
more scatter in the tubh — rnh relation than is observed. 
Fig . El shows t hat to be consistent w i th the tightest observa- 
tional limits Wemaine_ei 'al] l2002l: iMarconi fc Hun^ 120081 
iHaring fc R.ix ]20^4^ requires Upf < 450-650 km s“^. This 
constraint is consistent with the kick velocities pr e dicted 
by the general relativisti c analysis of iFavata et al.l (l2004^ 
an d iBlanchet et all (l2005h . Even so, if the prefactor velocity 
lies in the range Upf = 300 - 500kms“^, as these calculations 
suggest, then gravitational recoil provides a substantial con¬ 
tribution to the scatter in the ttibh — rnh relation. 

In models where discs are assumed to be stable to the 
formation of a bulge, all the observational estimates of the 
scatter lie well above the values that result from recoiling 
black hole ejections. In this case, we conclude that even kicks 
with velocities exceeding the range expected from current 
calculations do not produce a significant contribution to the 
scatter in the rriBH — rnh relation. There are two reasons for 
such a large difference between the two models. Firstly, when 
discs are assumed to be stable there are far fewer (~ 10%) 
bulge-bulge mergers and far more disc-bulge and disc-disc 


Figure 3. The scatter in the ruBH — mb relation as a function of 
the kick velocity prefactor, Upf. The curves show our estimates of 
the width of the distribution of MBn/^ideal BH induced by black 
hole ejections as discussed in the text. The solid and dashed curves 
represent ergs and age respectively. The thick lines show the esti¬ 
mates when disc instability produces bulges, while the thin lines 
represent models when this route for bulge formation is ignored. 
The horizontal lines show estimates of the scatter from various 
observational studies. Note that Tremaine et al. (2002) give only 
an upper limit, while Marconi & Hunt (2003) present two esti¬ 
mates to which they refer as conservative and unconservative. 

mergers involved in the formation of the galaxy. While the 
total number of mergers is the same in both models, the 
effect of the gravitational recoil depends only on the number 
of bulge-bulge mergers because, in our model, black holes 
only grow if there is a stellar spheroid. 

The second reason for the difference between the two 
cases is that when discs are assumed to be stable, most ejec¬ 
tions occur at early times {z ~ 2) when the bulges were 
smaller. The bulge then has enough time to regrow a black 
hole of the appropriate mass by subsequent galaxy mergers 
from which black hole ejection becomes increasingly diffi¬ 
cult. When unstable discs are assumed to form bulges, the 
majority of the ejections occur much later {z < 1) and sub¬ 
sequent mergers are unable to grow a large black hole. Note 
that no black holes are ejected at 2 < 0.7 and 2 < 0.3 for 
disc-stable and disc-unstable models respectively. 


5 THE BLACK HOLE DISTRIBUTION FROM 
N-BODY SIMULATIONS 

In this section, we use high resolution N-body simulations 
of galactic halos to track the location of subgalactic and 
galactic black holes as the halo grows by mergers. The sim¬ 
ulations provide spatial information that is not available in 
Monte-Carlo merger trees. We first describe the simulations 
and the way in which we model black hole escapees. 
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Ntot 

(10®) 

Nh, 

(10®) 

7?vir 

(/i-i Mpc) 

Nvir 

(10®) 

ghl 

14.6 

12.9 

0.110 

1.07 

gh2 

18.1 

16.2 

0.131 

1.74 

gh3 

18.0 

16.2 

0.170 

3.73 

gh6 

25.5 

22.2 

0.169 

3.76 

gh7 

19.2 

17.3 

0.156 

2.99 

ghlO 

13.4 

12.1 

0.133 

1.86 


Table 1. Parameters for the six A^-body halo simulations. The 
columns give: (1) halo label; (2) total number of particles in the 
simulation cube; (3) number of high resolution particles; (4) virial 
radius of the halo in h~^ Mpc defined as the distance from the 
centre to the radius at which the mean interior density is 200 
times the critical density; (5) number of particles within the virial 
radius. All halos were simulated in a cube of comoving length 
35.325 h~^Mpc in a ACDM universe, with a particle mass of 
2.64 X 10^h~^MQ in the ‘high resolution’ region. 


5.1 Galaxies in the N-body simulations 

We have used six Wbody simulations of galactic-size dark 
matter halos of final mass ~ 10^^ Mq to study the demo¬ 
graphics and spatial distribution of their black hole popu¬ 
lation, including black hole ejectio ns. The simulations w ere 
performed using the GADGET code iSoringel et alj[2QQlll in 
a flat ACDM universe with Dm = 0.3, h = 0.7, as = 0.9. 
These simulations have been previously studied by other au- 
thors JPowe^^^ ]|2nn3| [ Havashi et al ] l2004l : lNavarro et alJ 

120041 : 11] he,skind et a,]Jl2nnfJ^ and we refer the reader to these 

papers for technical details. Briefly, each simulation fol¬ 
lows the formation of structure in a cube of comoving side 
35.325 /i“^Mpc with a Lagrangian ‘high resolution” region 
around the halo of interest in which the particle mass is 
2.64 X h~^M q and ‘low resolution’ particles elsewhere. 
Tabled summarises the parameters of our six simulations. 

Halos and their su bstructures were ide ntihed using 
the algorithm SUBFIND JSpringel et al.l 1200 ill . The algo¬ 
rithm first identifies “friends-of-friends” groups using a link¬ 
ing length of 0.2 times the mean interparticle separation 
which approximately select s particles lying w ithin the viri- 
alized region of the halo llDa.vis et al.l ll 98, 5^ . Only halos 
with more than 10 particles, corresponding to a mass of 
2.64 X 10^h~^M q, are considered. Then, using an excur¬ 
sion set approach, SUBFIND identifies self-bound subgroups 
within each friends-of-friends halo. 

We use the method described by iHarker et al.l ll2005^ 
to construct full merger histories for the dark matter ha¬ 
los. Progenitor and descendant halos are identified at every 
timestep and tracked throughout the simulation in order to 
build the merger tree. The semi-analytic galaxy formation 
code GALFDRM (e.g see ICo1e et al1l2nnflt iRenson ^ 
iBaugh et alJl993i is then applied along each branch of each 
merger tree to obtain the properties of the central galaxy in 
each halo and its orbiting satellites. If the subhalo that hosts 
a satellite survives inside the parent halo, the position of the 
satellite is identified with the most bound subhalo particle. 
Some subhalos, however, are disrupted by tidal forces as 
they sink by dynamical friction and can no longer be iden¬ 
tified by SUBFIND. In this case, the satellite is placed at the 
centre of mass of the particles that made up the subhalo 
at the last time it was identified. If the harmonic radius of 
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Figure 4. The mean interior radial dark matter density profiles 
for six simulated halos plotted against radius in units of the virial 
radius, and normalised to the value of the mean interior density 
at the virial radius (long dashed curve). The solid curves show the 
mean interior number density of satellites as a function of radius, 
normalised to the value at the virial radius. The dotted curve 
shows the mean interior density of all black holes (both those 
residing in satellites as well as the wandering ones), normalised to 
the value at the virial radius, assuming a kick velocity, Upf = 300 
kms~^. 


these particles becomes greater than the distance between 
the satellite and the centre of the parent halo, the satellite 
is deemed to have merged into the central galaxy. 

Fig. 0] shows the spherically averaged dark matter den¬ 
sity profile of our six halos normalised to the mean dark 
matter density within the virial radius. The galactrocen- 
tric distance is plotted in units of the virial radius (see 
Table 1). We also plot the corresponding number den¬ 
sity profile of satellite galaxies with a V band magnitude 
Mv < —7. The dark matt er pr ofile follows the NFW 
fomjNavarro^^enk_^ WhitdllQQd Il997l quite closely (see 
iNavarro et al.l2004^ . The radial distribution of the satellites 
is shallower than that of the dark matter in the inner parts of 
the halo. This result is broadly in agreement with the density 
distributions of substructures previousl y obtained from high 
resolution N-body sim ulations (see e.g iGhinga et al1ll99Sl . 
l200(l iGao eIaill2004^ . Note, howeyer, that the profiles in 
Fig. El refer to satellite galaxies, not substructures and that 
19.3% of the satellites are not attached to a subhalo. Ney- 
ertheless, the radial profile of the satellites is similar to that 
of subhalos. 


5.2 The black hole population 

The growth of galactic bulges and their associated black 
holes is calculated along each branch of the merger tree 
as described in Section 3.3. At the final time, each central 
galaxy contains a black hole typically of mass 10^h~^M q for 
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Figure 5. The spatial distribution of mass and black holes in simulation gh2. The top two panels correspond to the model in which 
unstable discs produce bulges, while the lower two panels correspond to the model in which this bulge formation channel is ignored. The 
left-hand and right-hand panels correspond to kick velocity prefactors of Upf = 300 kms~^ and Upf = 500 kms“^ respectively. Black 
holes are represented by symbols whose area is proportional to their mass. The open circles denote black holes that lie at the centre of a 
satellite galaxy, while open squares denote wandering black holes that were ejected during a merger. The central galactic supermassive 
black hole is in the middle of each panel. The virial radius of the halo is indicated by the dashed circle. The subset of intergalactic black 
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models where unstable discs form bulges and Mq for 

models where disc instability is ignored. There are, in addi¬ 
tion, two other populations of black holes. Firstly, there are 
“galactic black holes” that reside in the bnlges of satellites 
in the halo; secondly, there are black holes that have been 
ejected from their host galaxy and which we term “wander¬ 
ing” black holes. The population of wandering black holes, 
in tnrn, is made np of “intrahalo black holes” that are still 
inside the host halo’s virial radius at 2 = 0 , and “intergalac- 
tic black holes” which have been ejected from the halo and 
lie outside the virial radius at 2 = 0 . 

5.2.1 Tracking black holes 

Black holes that remain within their host bnlge are always 
assumed to be located at the centre of the galaxy. When a 
black hole is ejected from its galaxy as a result of a recoil, we 
search for a particle in the simulation in the neighbourhood 
of the galaxy with appropriate velocity which, we assume, 
tracks the orbit of the black hole. To hnd as close a match as 
possible, we dehne a cost function and attach the escaping 
black hole to the particle that minimises the cost. The cost 
function we adopt is: 



where Ar^ and Vi^f. are the position and radial velocity of 
the ith particle relative to the initial black hole position and 
velocity, Vej = (Ukick “ is the velocity the black hole 

has when it escapes from the galaxy and r* = 0.33 Mpc. 

For modest kick velocities, Upf < 300 kms“^, a suit¬ 
able particle is usually found and 75% of the ejected black 
holes are identified with particles with Ci < 0.5. As Upf 
is increased, the identification becomes increasingly diffi¬ 
cult. Black holes with such large kick velocities, however, 
are likely to be ejected from the halo in any case. Thus, 
if no particle with Ci < 0.5 is found, we assume that the 
black hole has been lost from the halo to become part of a 
population of intergalactic black holes. 

The general properties (number density, luminosity 
function) and spatial distribution of the satellite popula¬ 
tion depend somewhat on whether or not we assume that 
unstable discs generate galactic bulges. The statistics of the 
associated black holes, on the other hand, depend strongly 
on this assumption since many more bulges and black holes 
are formed when disc instability is taken into account. As 
above, in what follows, we take the unstable disc case as 
our fiducial model, but summarise also results when only 
mergers are assumed to give rise to bulges. 

5.2.2 The distribution of black holes 

The spatial distribution of black holes and dark matter in 
one of our simulations is illustrated in Fig. The top two 
panels correspond to the model with unstable discs and the 
lower two to the model without unstable discs. In both cases, 
the left-hand column shows results for a kick velocity prefac¬ 
tor of Upf = 300 kms“^ and the right-hand column for a kick 
velocity prefactor of Vp( = 500 kms“^. The greater efficiency 
of bulge formation in the unstable disc case is reflected in 
the larger number of black holes in this model. Those black 


holes that are still associated with a satellite galaxy are in¬ 
dicated by circles while those that have been ejected (the 
“wandering” black holes) are indicate by squares. 

Tests of the mass resolution of our calculation in which 
we artfficially increased the minimum halo mass used in 
the merger trees indicate that our catalogues are essen¬ 
tially complete for satellite galaxies with V-band luminos¬ 
ity brighter than -7 (corresponding to a mass larger than 
1.4 X lO®/i“^M 0 ). To this limit we find that, for the unsta¬ 
ble disc model, each simulated galaxy halo at 2 = 0 contains, 
on average, 52 satellites within its virial radius. The median 
bulge stellar mass of this population is 6.5 x 10^ h~^M q with 
a corresponding black hole mass of 6.5 x 10^h~^M q in the 
ideal meH — rnh relation. The fraction of these satellites 
that retain a central black hole at the final time depends 
weakly on the assumed prefactor for the kick velocities. For 
Wpf $ 300 kms“^ 85.4% of satellites retain a central black 
hole at the final time, but this fraction is reduced to 84.2% 
for r’pf S 500 kms“^. 

The number density profiles of the black hole popu¬ 
lations of each of our simulations, assuming our fiducial 
Wpf = 300 kms“^, are shown by the dotted lines in Fig. |1| 
These can be compared with the dark matter density pro- 
hles and the corresponding number density prohles of satel¬ 
lite galaxies with a V band magnitude Mv < — 7. In order 
to compare wandering and satellite black hole populations 
of similar masses, for this plot we selected only black holes 
with mass greater than that of the black hole in the low¬ 
est mass bulge of the satellite sample. The density proHle 
of black holes is seen to be intermediate between that of 
the dark matter and that of the satellites and, in a cou¬ 
ple of cases, wandering black holes are found as far in as 
rvir/30. Our understanding of this behaviour is that these 
black holes are ejected from the progenitors of both present 
day satellite galaxies and the progenitors of galaxies which 
by the present have merged with the central galaxy. As the 
ejection velocities are not large, these wandering black holes 
initially have orbits similar to the galaxies that ejected them 
but, being lighter, they are not subject to dynamical friction 
and so, unlike the satellite galaxies, those near the centre are 
not removed by merging with the central galaxy. 

5.2.3 The blaek hole mass function 

In the case of the fiducial model (in which unstable discs 
form bulges), each Milky-Way mass halo generated on av¬ 
erage 80 black holes of mass greater than for 

300 < Upf < 500 kms“^ during the course of its formation. 
Of these, at 2 = 0 approximately half were wandering black 
holes and the other half were galactic black holes retained 
by their host satellite galaxies. Of the wandering black holes, 
the fraction that are ejected with sufficient velocity to com¬ 
pletely escape the halo depends more sensitively on Upf. For 
Upf = 300kms“^, 42% become intergalactic wandering black 
holes while for Upf = 500kms“^ this fraction increases to 
60%. These numbers and the corresponding ones for the 
model in which disc instability is ignored are given in Ta¬ 
ble 2 . 

At the present day, there are only 50 black holes more 
massive than 10^h~^M q within the central 100 /i“^kpc of 
the halo. Black hole passages through the disc are therefore 
rare and unlikely to have affected the structure of the disc 
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Table 2. The mean number of satellite galaxies and black holes for both the fiducial model in which unstable discs form bulges and 
the model in which disc instability is ignored. In all cases we only consider satellites within rvir with Mv < —7 and black holes with 
-^^BH > 10^h~^MQ. The first section of the table lists the mean number of satellites with bulges and their median and minimum total 
stellar masses. The second section lists, for both Upf = 300 kms“^ and Upf = 500 kms“^, the mean total number of black holes within 
the virial radius, Nbh^ of which NgatBH galactic black holes residing in satellite galaxies, and NjntrahaloBH wandering black holes 
within the main galaxy halo. The last row gives NintergalacticBH? mean number of wandering black holes that have escaped beyond 
the virial radius of the main halo. 
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Figure 6. The total mass functions for halo black holes. The 
left hand panel shows the average number of black holes per halo 
for Upf = 300 kms“^ while the right hand panel shows the same 
quantity for Upf = 500 kms“^. The upper two panels show his¬ 
tograms for models where discs are unstable to bulge formation 
while the bottom two panels show data for models where insta¬ 
bility is ignored. The red dashed histogram represents the mass 
function of black holes that reside in satellite galaxies; the blue 
solid histrogram represents the mass function of wandering black 
holes that have been ejected from their host bulges. The black 
histogram is the combined mass function. 


substantially. A large population of intergalactic black holes 
exists beyond the virial radius of the halo. Some of these 
are attached to small galaxies that will eventually become 
satellites, but others are black holes that were ejected from 
the halo altogether. (Only the fraction of these that were 
tagged with a particle are shown in Fig. [3) 

Fig. 1^ shows the mass function of halo black holes, av- 


Figure 7. The fraction of black holes that remain in satellites as 
a function of black hole mass. The left two panels are for Upf = 
300 kms“^, the right two are for Upf = 500 hms"*^. The upper 
two panels are for unstable disc models, while the lower two panels 
are for stable disc models. 


eraged over all six simulations, for two assumed values of 
Wpf in the models with and without unstable discs. Black 
holes in satellites within the virial radius whose host galaxy 
is brighter then My < —7 and all wandering black holes 
are included. In all cases, the central supermassive black 
hole of mass ~ 10® — 10^ Mq has been omitted from the 
sample. In the unstable disc case, the medians of the dis¬ 
tribution are and 10^ ®®/i“^Mq for Upf = 300 

and 500 hms”*^ respectively. In both cases, we expect a typ¬ 
ical halo to have ^ 3 black holes of mass M > 10®/i“^Mq. 
The mass distribution is slightly wider for the larger value 
of Upf, reflecting the larger efficiency of black hole ejection 
in this case. In the case where disc instabilities are ignored, 
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their are fewer black holes (note the difference in the scales 
of the y-axes in the upper and lower panels of this figure), 
but the distributions are shifted to slightly higher black hole 
masses with a median of ^ Mq and virtually no 

black holes of mass M > Mq. The reason why the 

median is shifted to slightly higher masses in the stable disc 
case is because in this case there are fewer ejections, thus 
allowing satellite galaxy black holes to grow larger. 

Also shown in Fig.|^is the result of splitting the black 
hole mass function into the two populations: the satellite 
galaxy black holes and wandering black holes (red dashed 
and blue solid line respectively). The bimodal nature of the 
total mass function is visible as the result of these two dis¬ 
tinct populations. The median of the satellite black hole 
mass distribution is ~ Mq while the wandering 

black hole population peaks at ~ Mq for unsta¬ 

ble disc modes. In the stable disc models, the bimodality 
is also visible although slightly less pronounced. This ten¬ 
dency for the wandering black holes to be more massive than 
those retained in the satellites is shown clearly in Fig. |71 
which plots the fraction of black holes of each mass that 
are associated with satellite galaxies as a function of black 
hole mass. We see that below 10® Mq the fraction associated 
with satellite galaxies is a decreasing function of black hole 
mass. For low masses (M < h~^M q), more than 90% 

of black holes are in satellites. However, at higher masses 
(10®'® < M < lO®/i“®M 0 ), the majority of black holes are 
found to be wandering. The reason for this bias in the masses 
of wandering black holes compared to satellite black holes is 
the fact that larger satellites are more likely to merge and 
experience a kick then smaller ones. Since the dynamical 
friction ti mescale is inverse ly dependent on satellite mass 
(eq. 4.16 in ICole et al.ll200(ill . we expect the physics of kicks 
to have a larger effect on the higher mass satellites and 
bulges than on the lower mass ones. The lower mass satellites 
will have experienced few if any ejections and their central 
black holes will be unaffected by velocity kicks. 

In principle, a large population of wandering black holes 
could affect the match between the total mass in black holes 
at the present day and the m ass in quasar remnants inferred 
from energy considerations ^SoltarJll982^ . In our calcula¬ 
tions, however, wandering black holes make up, on average, 
only 2.6% and 3.9% of the total black hole mass in our sim¬ 
ulated galactic halos for Wpf = 300kms“® and 500kms“® 
respectively. Wandering black holes do, however, make up a 
large fraction of the intrahalo black hole population (ie. all 
halo black holes excluding the central one). For these two 
values of Upf, unattached wandering intrahalo black holes 
make up 31.8% and 48.4% of the total intrahalo satellite 
black hole mass. 

Another noteworthy feature of the black hole mass func¬ 
tions of Fig.l^is the small population of very massive black 
holes with Mbh > 10®/i“®Mq. This population is primar¬ 
ily composed of wandering black holes, although a portion 
of it in the unstable disc model consists of black holes in 
massive satellites whose bulge masses are ~ 10^'^'^h~^M q, 
as evidenced by the high mass peak in the upper panels of 
Fig.lZl In stable disc models, this population is composed 
entirely of wandering black holes. These massive wandering 
black holes originate from ejections from the central galaxy 
that occured when the central galaxies bulge had a mass 
of ~ 10^h~^M q. The central galaxy has since undergone 


subsequent phases of bulge formation to bring its mass to 
~ and, in the process, has regrown a central 

black hole of mass < h~^M q. 


6 DISCUSSION AND CONCLUSIONS 

The gravitational recoil of merging black holes is an im¬ 
portant physical effect in a universe built up hierachically 
through the repeated merging of galactic units. Whenever 
galaxies that host black holes merge, the black holes them¬ 
selves will coalesce and there exists the potential for the 
remnant black hole to be ejected from the galaxy, provided 
the recoil velocity is high enough and the galactic potential 
shallow enough. 

We have examined the role that gravitational recoil 
plays on the demographics of black holes. We find that the 
process of ejecting black holes from galaxies is efficient if 
bulges and black holes grow both in galaxy mergers and 
as a result of discs becoming unstable, because in this case 
nearly all bulges contain a black hole. In models where disc 
instability is ignored, this process is not as efficient because 
fewer bulges and associated black holes exist and because 
black hole - black hole mergers tend to occur late when the 
galactic potential wells are deeper. In the former case, black 
hole ejections produce a significant contribution to the scat¬ 
ter in the mBu — rrih relation. In fact, conservative estimates 
of the scatter in the observed relation constrain the recoil 
prefactor velocity to Vp{ < 500 kms“®, which is consistent 
with the g eneral relativiti s tic ca lculations of iFavata et alJ 
^2004^ and iBlanchet etld] ||200^. 

Is there any empirical evide nce for the kind of blac k hole 
processes present in our model? ICoccato et al.l (l 200 fjl have 
recently reported the discovery of a black hole in NGC 4435 
whose mass is smaller than about 20 % of the value expected 
from the mBH — rrih relation. Black holes with a smaller 
than expected mass arise naturally in our model, although 
masses as extreme as this could be rare. For example, in our 
model with Upf = 300 kms“®, only 2 — 3% of black holes 
have a mass that deviates as much or more from the mean 
relation as that of the black hole in NGC 4435. However, if 
Upf = 500 kms“® this fraction raises to ~ 20%. 

Indirect evidence for black holes with masses below 
those exp ected from the itibh — n ih rel ation has been pre¬ 
sented by IColbert fc Mushotzkvl lll999ll . They argue that 
the compact X-ray sources often seen near the centre of 
elliptical and spiral galaxies could be black holes of mass 
~ 10^ — 1 O"^M 0 . Interestingly, these sources are often dis- 
placed from the centre by a few hundred parsecs. Similarly, 
iNeff et alJ J2003h have discovered a group of off-centre com¬ 
pact X-ray sources in the merger remnant NGC 3256 which, 
they argue, could be intermediate-mass black holes. In our 
model, these objects might be identified with the black holes 
of infalling satellites or with recently merged black holes that 
have been kicked out of the galactic centre, their growth 
stunted as a result. 

To investigate the spatial distribution of black holes in 
Milky-Way like galactic halos we used a set of N-body sim¬ 
ulations to track both the satellite galaxies that host black 
holes and also the black holes that are ejected from their 
host galaxies. We find that the black hole mass function is 
bimodal, being composed of two overlapping populations. 
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The lower mass population consists of black holes at the 
centre of orbiting low mass satellites that have not under¬ 
gone recent mergers, while the slightly higher mass popula¬ 
tion is composed of wandering black holes that have been 
ejected from mergers that formed the central galaxies and 
the more massive satellites. Among the latter population, 
we find a few supermassive (< 10®h“^M©) black holes that 
were ejected from the main progenitor of the central galaxy 
sufficiently early such that the bulge of the central galaxy 
has had enough time to regrow a sizable black hole. There 
is also an intergalactic population consisting of black holes 
whose recoil velocity was large enough not only to unbind 
them from their host galaxy but also from its halo. 

In the future it may be possible to detect the formation 
of wandering black holes directly using instruments such as 
LISA^ that can directly measure the gravitational radiation 
emitted in the black hole - black hole merger that ejects 
the remnant. In the meantime, detecting these black holes 
presents an interesting challenge. In practice, an ejected in¬ 
tergalactic or intrahalo black hole is likely to bring along 
a small cusp of stars tightly bound to it as it escapes the 
halo’s potential. These stars would provide the only mea¬ 
surable way of detecting such black holes. Most likely, there 
would not be any gas that could be accreted and hence the 
black holes would not be observable in the conventional way. 
However, as the orbits of the stars which the black hole 
brought with it from the bulge decay, the stars may plunge 
into the black hole. In addition to emitting a burst of grav¬ 
itational radiation, this type of infall would tidally disrupt 
the star and create a small accretion disc which would radi¬ 
ate according to the standard physics of accretion discs. The 
gravitational radiation signature and accretion disc emission 
wou ld provide the best w ay to identify these black holes. 

iMagain et al.l ll200?jl have recently discovered a quasar, 
HE0450-2958, which has no visible host galaxy and, at 
first sight, is a good candidate for an escaped black hole. 
iHaehnelt. Davies fc Ree^ ^200!^^ have indeed suggested that 
the quasar may have been ejected from a nearby ultralumi- 
nous infrared galaxy (ULIRG) either through the gravita¬ 
tional recoil process we are considering here or through a 
gravitational slingshot associated with three or more black 
holes inv olved in the merg e r resp onsible for the ULIRG. 
However, iHoffman fc LoebI feOORIl have argued that the 
quasar is much too far from the companion galaxy to have 
been ejected with a velocity of 300 kms“^ and so favour a 
slingshot mechanism. 

Ultraluminous X-ray sources (ULXS) or micro-quasars 
have long been regarded as candidates for intermediate- 
mass black holes, Mbh ~ 100 — lOOOM©, r adiating near the 
Eddington limit (see e.g. iFabbian^llOSOt iMushotzkvll2004l 
and references therein). This is exactly the kind of object 
that would be naturally identified with the orbiting intra¬ 
halo black holes predicted by our model. This interpretation, 
however, has two difficulties. Firstly, ULXs tend to be asso¬ 
ciated with star-forming regions and their frequency seems 
to be correlated with the galactic star formation rate (see 
e.g. IWardlliooil : iKaaret et al ] l2004l and references therein). 
These facts lend support to the view that ULXs are stel¬ 
lar mass black holes emitting beamed radiation at highly 


http://lisa.jpl.nasa.gov/ 


super-Eddington rates. The second arguement against iden¬ 
tifying ULXs with the intrahalo black holes in our model is 
the difficulty of finding a suitabie source of materiai for the 
biack hoie to accrete. There are, howev er, some exampfes o f 
ULX’s, most notably the ULX in M82 iKaaret et al.ll20oili . 
that are probably much too bright to be explained even as 
exotic super-Eddi n gton luminosity s tellar mass black holes 
iKing et alJl20oill . iKing fc DehnenI i200fill call objects like 
this “Hyperluminous X-ray sources” (HLX) and argue that 
these objects are precisely the intrahalo black holes associ¬ 
ated with satellites in our model which switch on when they 
come close to the galactic centre. Since the exact mechanism 
by which these black holes would be activated is uncertain, 
we cannot predict how common this phenomenon should be. 
However, our model contains, in principle, a plentiful sup¬ 
ply of intermediate mass black holes orbiting in the halo of 
every galaxy which is sufficient to account for the presence 
of a few ULXs or HLXs in most gaiaxies. 
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